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Abstract 

By the weight of a Boolean function /, denoted by wt{f), we mean the number of inputs 
for which / outputs 1. Given a promise that an n- variable Boolean function (available in the 
form of a black box and the output is available in constant time once the input is supplied) is 
of weight either wN or (1 — w)N {0 < w < 1, N — 2"), we present a detailed study of quantum 
algorithms to find out which one actually it is. To solve this problem we apply the Grover's 
operator. 

First we consider the restricted problem. Given a promise that an n- variable Boolean func- 
tion is of weight either [iVsin^ IFPlfl L-^cos^ SFPrfl (L^l means the nearest integer cor- 
responding to the real value q), we show that one can suitably apply Grover's operator for 
/c-many iterations to decide which case this is with a probability almost unity for large n and k 
in 0{poly{n)). On the other hand, the best known probabilistic classical algorithm has a suc- 
cess probability close to 0.5 (from above) after k many steps when k is large. We further show 
that the best known probabilistic classical algorithm can achieve a success probability almost 
unity only after A:" many iterations where s > 2. This indicates a quadratic speed up (and 
also agrees to the quadratic speed up by the use of Grover's algorithm in database search) on 
time complexity in the quantum domain with respect to the best known result in the classical 
domain. 

Second, we modify the basic randomized algorithm into a sure success algorithm, which can 
distinguish Boolean functions of weights wN or (1 — w)N for any w, {0 < w < 1). To do that we 
have exploited a sure success Grover search algorithm, which modifies the very last operation. 
For the weight decision problem, we show that the very last two operations should be changed to 
distinguish any weight with certainty and found the phase conditions for the last two operations. 

As quantum counting methods exist, which can count the number of solutions, here we 
compare our method with that. Since the quantum counting method needs to exploit period 
information, which requires many Grover operations, we have found that our method is faster 
than the quantum counting method. 

Keywords: Quantum Algorithm, Boolean Function, Grover's Operator, Weight Decision Prob- 
lem. 
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1 Introduction 



David Deutsch designed a quantum algorithm, which evaluates whether the two outputs of a 
Boolean function is the same or not using only one function evaluation . Deutsch- Jozsa gen- 
eralized the Deutsch's algorithm for more general case such as whether the Boolean function is 
constant or balanced "W. Deutsch- Jozsa algorithm had been proposed to show an exponential 
speed-up in quantum machine than the classical machines. The most important contribution in 
this area has been achieved when Peter Shor discovered a polynomial-time quantum algorithm 
for factoring and computing discrete logarithms, which is exponentially faster than the classical 
best known algorithm jl7j . Because this quantum factoring algorithm works in polynomial time 
compared to exponential time in classical method, many researchers started to find other applica- 
tions. Lov Grover discovered a quantum database search algorithm, which is quadratic faster than 
classical database search algorithm Meanwhile, database search algorithm is one of the most 
widely used algorithm in the computer applications, the impact is so huge and more researchers are 
interested in other applications of quantum database search algorithm. This research also focuses 
on the applications of Grover database search algorithm. 

In this research, we assume that there is a promise in the Boolean function as it has a weight wN 
or (1 — ■w)N. Initially, for better understanding, we formulate the problem for some special weight 
cases such as [N sin^{ 2j^^i §)] or [N €05^( 2^^;^ f)] , by applying k many Grover operations. Here the 
quantum algorithm presented is of randomized in nature, but the success probability is arbitrarily 
close to unity. Next we consider the general case to distinguish functions having weight either wN 
or (1 — w)N and we also consider a deterministic algorithm. This requires changes in the very 
last two Grover operations with phase conditions. Briefly, from l*** to {k — 2)*'* steps, the original 
Grover operators are used, but, in the last two steps, two different Grover operators are used with 
the phase conditions. Also we found that the phase conditions depending on the required number 
of Grover operations. Meanwhile, because quantum counting algorithm was already proposed, we 
compared two methods. Since the quantum counting method requires more Grover iterations to 
find period information, we can conclude that our method is faster than this method. 

2 Preliminaries 

A Boolean function f{xi, . . . ,Xn) on n variables may be viewed as a mapping from {0,1}" into 
{0, 1}. A Boolean function / is constant if f{x) = c for all x G {0, l}", c G {0, 1}. That means 
wt{f) is either or A^. A Boolean function / is balanced if wt{f) = N/2. 

Given a promise that the function / is either constant or balanced, one may ask for an algorithm, 
that can exactly answer which case it is. Note that throughout this document we consider that 
any Boolean function / is available in the form of an oracle (black box) only, where one can apply 
an input to the black box to get the output. A classical algorithm needs to check the function 
for A^/2 -|- 1 inputs in the worst case to decide whether the function is constant or balanced. It is 
known that given a classical circuit for /, there is a quantum circuit of comparable efficiency which 
performs a transformation Uf that takes input like \x,y) and produces the output \x,y(B fix))- 
Given such a Uf, Deutsch- Jozsa ;6. provided a quantum algorithm that can solve this problem in 
constant time, indeed, in a single evaluation of Uf. The circuit for their algorithm is given in Fig.^ 

Algorithm 1 Deutsch-Jozsa Algorithm jB] 
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Figure 1: Quantum circuit for Deutsch-Jozsa Algorithm 
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Measurement at M: all-zero state implies that the 
function is constant, otherwise it is balanced. 



The Deutsch-Jozsa algorithm yields an exponential speed-up relative to any exact classical 
computation. This provides a relativized separation between EQP and P with respect to the oracle 
/ (see |16| for basic notions of complexity theory). 

Now we discuss a constant time quantum algorithm to distinguish Boolean functions of weight 

N/4 and 3iV/4 [7|. We replace iJ®" by Grover's matrix at the output side of Uj in Fig. Qto get 

a circuit shown in Fig. |2l and show that this solves the problem. In 2001, Green and Pruim [Jj 

NP 

presented a relativized separation between BQP and P using a nice technique based on Grover's 
algorithm 0. Green and Pruim's work relied on a complexity theoretic formulation, whereas our 
analysis here is directly related to weights of Boolean functions. Note that a similar question has 
been discussed in [IJ Section 5]. There also the problem was not exactly posed as a discrimination 
problem, but as a search problem. 

We denote the iV x Grover's matrix G„ as G„ = //®"(2|0)(0| - 1)/?®" = ^ Y.x,y " 1- 

It is known that this operation may be constructed with 0(log A^) quantum gates jl^. The circuit 
is shown in Fig. [2 and the steps of the algorithm are as follows. 



Algorithm 2 Randomized Algorithm to distinguish wt{f) = 2>N/A or N/A 
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Figure 2: Quantum circuit to distinguish wt{f) = N/A and wt{f) = SN/i 
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4. m = GM)^ 

5. Measure the resulting state Gn\tp'2) the 
computational basis and let the result be x. 

6. if f{x) = then wt{f) = 3iV/4 else wt{f) = N/A. 



The key result proving the correctness of this algorithm is as follows and we present a proof as it 
will be discussed in details in the following section. 

Theorem 1 = Ex:/(x)=o ^^^i^W^\^) + Y.x:f{x)=i and Algorithm H produces 

a correct result. 



EJ^)-E.(-lK^")k>)- From 



Proof: GM',) = E.e{o,i}"(-l)'^^'^ 1^) = (^^^^^ 

which the result follows. 

If wt{f) = N/A then the probability amplitude of all the |x) for which f{x) = vanishes. So 
on measurement we will get some x for which /(x) is 1. On the other hand, if wt{f) = 3N/A then 
the probability amplitude of all the \x) for which f{x) = 1 vanishes. So on measurement we will 
get some x for which f{x) is 0. ■ 



3 Repeated Application of G rover's Operator 

Note that Grover's search algorithm ,Hj uses repeated applications of G„. Motivated by the same 
idea we now analyze in detail the repeated application of GrJJf. One may refer the important 
papers like|2llHllliniIiniIIllIl2lIiniIIl[II3IIHl where the Grover's operator has been used. We also 
make an elaborate study to present a comprehensive understanding of this problem in this section. 
We start with a modification of Algorithm [21 
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Algorithm 3 Randomized Algorithm for Weight Decision Problem 2>N/A and N/A 
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a:G{0,l}" ^ 



lV'o) = |0)®"|l) 
Let us denote iV^'i 

l^2) = f//(|^i)^) 

Let us denote first n qubits of \1IJ2) as |V'2)- 

m = G„|V^)^°^. Let W^) = GnW^). 
i = 2 + 1. 

If {i < k) denote iV'a) by 1-01) and go to step 4. 

Measure the resulting state ™ the 

computational basis and let the result be x. 

if k is odd and if f{x) = then wt{f) = [A^cos^ 

else ^^i(/) = [iVsin^^^fl. 

if k is even and if f{x) = 1 then wt(f ) = [A^cos^ 

else wt{f) = [iVsin^^l]. 



2fc+l 2 



2fc+l 2 



In this section we will consider a number of iterations k > 1 and then show how we can distinguish 
the weights with a very high (almost unity) success probability. Without loss of generality we detail 
the analysis by taking odd k. 



Theorem 2 Let denotes the amplitude of the states \x) where f{x) = and denotes the 
amplitude of the state \x) where f{x) = 1 after k-th iteration. Then, = {2 ^~'^^^'^ — l)afc_i — 



' N 



Jk-1 



, h 



2^^a,^, + (2 



N-wt{f) 
N 



l)bk-i, with initial conditions gq = bo = 



Proof: The proof follows from Algorithm |21 in a similar fashion as the proof of Theorem ^ ■ 
Our interest is to investigate the zeros of ak and b^. It may be noted that the solutions to the 
recurrence relations are given by 



1 sin(2A; + l)/3 
O'k = : — jOfc 



1 cos(2fc + l)/3 
cos (5 



where sin^ fi 



N 



sin/3 



(= u, say for notational convenience). Note that this recursion and Theorem[21 



have been described in Clearly, the factor does not play any part in determining the 

zeros of and bk. The zeros of ak and bk are given by 



sin(2A; + l)(3 = 0, sin /? / and cos(2A; + 1)/? : 

respectively Now sin(2A; + l)/3 = ^ (2/c + l)/3 = /vr ^ /3 

C 

(2fc+l) 2 



0, cos /? / 

where / G Z. Also 



(2fe+l) 

cos (2/c + l)/3 = ^ (2A; + 1)/? = (2m - l)f ^ /3 : 

As we are interested in distinct roots u = svo? (3 of = (respectively 6^ = 0), it is clear that 
we will get the distinct roots when \ < I < k (respectively 1 < m < fc). We can summarize the 
above discussion in the following result. 
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Proposition 3 The k distinct roots of ak = and bk = are sin^ (2k+i) '^^^ ^(2k+i) f i^^spec- 
tively where 1 < l,m < k. 

Proposition 4 For each root of the equation = there is a corresponding root of the equation 
bk = so that their sum is 1. 

Proof: The roots of ajt = and 6fe = are of the forms siii^ (2fc+i) '^^^ s\v? ^(^2k+i) 5 ^^spectively 
where 1 < l,m < k. Let us consider the pairings of /, m such that I + m = k + 1. Now sin^ [2k+i) ~^ 

sin^ ^^^Tfi ~ '^hich gives the proof. ■ 
Given any k, Algorithm |31 can distinguish whether a Boolean function is either from the weights 
[N sin^ ^{2k+i) f "I from the weights [A^ cos^ ^{2k+i) f "I ' '^i^h good success probability for m = 
1, . . . ,k. However, here we are interested in distinguishing two Boolean functions which are closest 
in weight to balanced functions, i.e., of weight ^ = 2"~^. 

Definition 1 Let ra^k be the root of ak = such that \ra^k ~ 0.5| < \pa,k ~ 0.5| for any root Pa^k of 
ak = 0. Similarly let r^^k be the root of bk = such that jrfe^fc — 0.5| < \ph^k — 0.5| for any root ph^k 
of bk = 0. Let us denote pk = niinlr^ fc, r^^fc}. 



Proposition 5 

Pk 



. 2 2a — 1 vr 

r„ k = sm , when k = 2a — \ and m = a 

4a - 1 2 ' 

2 2a vr 

= Thk = sin , when k = 2a and I = a. 

'^''^ 4a + 1 2 ' 

Proof: Let A; = 2a — 1. For m = a, the root of = is pa.fc = sin^ la-i f • ^'-'^ m = a + 1, 
the root of = is p'^ ,. = sin^ . As 2a + 1 > 2a, sin^ > sin^ , which gives, 

■ 2 2a+l TT I 2 2a TT ^ 1 ■ - 2 2a+l tt , - 2 2a-l n ^ -i • „i„2 2a+l tt _ 1 ^ 1 _ ■ 2 2o^ tt 

4a~l 2 ^^^^ 4a-l 2 ^ ^-^'^ 4a-l 2 ^''"^ 4a~l 2 ^ ^' l"^'' ''"^ 4a-l 2 2^2 ''I" 4a-l 2 

which gives ^ — ^ > ^ — /Oa,A:- Since all the other roots of Ofc = are either less than pa^k or 
greater than p^^, we get r^ fc = pa^k- It is also clear that rh^k = 1 — ^a,k and hence here r^^fc > ^a,k- 

So /Xfc = ra^k- 

Let A; = 2a. For m = a, the root of 6^ = is ph^k = sin^ Aa+i f • m = a + 1, the 

root of 6fc = is ;t = sin2 . As 2a + 2 > 2a + 1, sin^ f^f > sin^ f^f , which gives, 

■ 2 2a+2 TT _i_„„„2 2a+l TT ^ i ■ - 2 2o+2 tt _i_„;„2 2a tt ^ i • - 2 2a+2 tt _ 1 ^ 1 _h;t,2 2a tt 
''"^ 4a+l 2 '^^^^ 4a+l 2 ^ ^-^'^ 4a+l 2 "^^"^ 4a+l 2 ^ ^' l"^'' ''"^ 4a+l 2 2^2 ''"^ 4a+l 2 

which gives P;, fc — | > ^ — P6,fc- Since all the other roots of 6^ = are either less than pb^k or greater 
than p'^ we get r^^k = Pb,k- It is also clear that r^^fc = 1 — r^^fc and hence here ra^k > ^b,k- So 
/^fc = »^fe,fc- ■ 

Theorem 6 pk = sin^ ^fc+i f '^'^^ /^fc < /^fc+i < 0.5. 

Proof: From Proposition [3 it is clear that pk = sin^ ^fc+i f • /^fc+i = sin^ 2 fcX 3 f • ^ow — 
IFTT = (2fc+3K2fc+i) > 0' which gives pk+i > Pk- Further it is easy to see that pk+i = sin^ < 
sin^ f < ^- ■ 
As ak and bk can be seen as polynomials in u, we now refer them as ak{u) and bk{u), respectively. 
It is clear that akipk) = 0. Now using Algorithm El we can distinguish two Boolean functions of 
weight pkN and (1 — Pk)^ ■ Unfortunately, pkN may not be an integer and in that case we have 
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to consider a Boolean function of weight /u'^A^, where ^'^^N is an integer and l^'^A^ — fikN\ < 0.5. 
Thus we win be using the Algorithm |S1 to distinguish between Boolean functions of weight fi'f^N 
and (1 — lJ''f.)N. This will incorporate some error in the decision process. However, we will show 
that this error is almost zero for large A^. 

Theorem 7 Consider Boolean functions on n variables and let N = 2^. After k iterations, k in 
0{poly{n)), the quantum algorithm (Algorithmic can distinguish two Boolean functions of weights 

[N siiP 2k+i f 1 ^"■^ [A^ cos^ 2k+i f 1 ^^^^ success probability > 1 — ^^^^''''^ which is almost unity 
for large N . 

Proof: We have = l^n /j^"*^ ~ ^' when sin^ /3 = Let //'^ = sin^/3'. We like to calculate 

the value of ak{fi[) = ^^^^^^Sj^- 

As \fi'j^N-fikN\ < 0.5, weget |sin2/3'-sin2/3| < ^. Thus (sin /?'+sin/3)| sin /5'-sin/3| < ^, i.e., 
(2sin^cos^)(2cosM|sin^|) < i.e., (2 sin ^ cos M)(2 cos ^| sin ^|) < 
^, i.e., sin (/?' + f3)\ sin {/3' - p)\ < ^. This implies, | sin - /3)\ < 2N sm\f3'+(3) - ^o^^ that 
i = ^1 < < i for fc > 1. Now = sin^p. So, f < /3 < f . As p ^ f3' , f3 + p' ^ 2/?. 
Due to the small difference between /3 and it may happen that on the lower side /3 + /3' may 
marginally be less than 2^ and at the higher side may marginally exceed 2j. Thus it is safe to 
assume sin (/3 + /?') > i. Hence | sin {P' - P)\ < ^ = jf- Since |^| < | sin (/?' - P)\ < \P' - P\, 

we can write | ^ 2 ^ I < Jj^ ~ P\ < W- 

One can take (l){P) = ^™ g^n^'*'^^ and use Taylor's series expansion for (f>{P) to get the upper 
bound on | i^n^ - ^^^^^K^^l- We consider (/>(/? + h) = ^^^^^^^r^, where h = P' - P is a 
small quantity. Now + /i) — (/)(/?) « /i '^^^^^ with error term bounded by Ri{P), the remainder 
when only the first term in the Taylor's series is considered. As k increases, the value of P falls 
in the neighbourhood of ^. It can be checked that | \p=^ < V2{2k + 2). Also we calculate 

1^1(1 )l = ¥l^l/3=f < ¥V2((2A: + 1)2-1)2 + 16, where < a < 1. As \cp{p + h)-m\ < 

I'i't'Wl I I P r^-ll rrt I sii^ (2fc+l)/3' sin(2fc+l)/3 | ^ 8(fc+l) 

|-3;a-|/3=f + we get | — ^ — | < 

Since ^'^^(2^+^)^ - we have ^^'^(2^+^^' < 8(^+1) Thus aJn' ) - 1 sm(2fc+l)/3' 8(fc+l) 
^^'^^i sin/3 - u, we nave ^j^^, < ^ . inus a^^/x^j - ^ ^j^^, < • 

Thus, if wt{f) = Hk^i then we would have got x such that f{x) = 1 with certainty. As fJ^tAf, 

may not be an integer, we have considered functions with wt{f) = ^'^N which is an integer such 

that \^'j^N — fikN\ < 0.5. In this case the probability of (wrongly) observing an x such that f{x) = 

is (a.(^l,))2(l - ^,',)N < if^m - l^dN < 

Similarly, if wt{f) = (1 — fJ-k)^^ then we would have got x such that f{x) = with certainty. As 
{1 — fik)Al , may not be an integer, we have to consider functions with wt{f) = (1 — /i^)A^ which is an 

integer. In this case the probability of (wrongly) observing an x such that f{x) = 1 is < ^^^^2^^ ■ 
■ 

Since the function / is available in the form of an oracle, the best known classical probabilistic 
algorithm can work as follows. For k many iterations it can present random inputs to the oracle 
and guess the function is of lower weight if the output zero appears more frequently and guess the 
function is of higher weight if the output one appears more frequently. As we consider the majority 
rule, we choose the number of iterations as odd in the classical probabilistic algorithm. This will 
always guarantee majority of either output zero or output one (in the case of an even number of 
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iterations there may be the possibihty of a tie and then a random decision has to be taken). For 
general analysis, the estimate of probabilities will remain almost the same, and we only present the 
analysis when the number of iterations is odd. 

Note that for classical lower bound proofs to determine the majority one may refer to 
However, we provide a complete result with proof which is suitable for our analysis. Here, the 
probability of the correct answer in a single step is cos^ 2fc+i f • After g{k) many iterations, g{k) in 

n{k), probability of success Ps = I - E»3^ (^^f )(cos2 ^^^^Y{sm'^ Let us now 

present the following technical result. 



Proposition 8 For odd positive integer k and g{k) be in Q{k), let 



g(fc)-i 

2 



E{k,g{k))= (tO(-^'^7^rTTT)^(^i-'^^^'^'^"^ 



i=0 

Then 



2{2k + l)' ' 2(2A; + 1) 



limk^oo E{k,k) = 0.5, 
Ihnk^ oo E{k,k'^) > 0.2 and 
limfc^oo = for s > 2. 

Proof: Let X/^. follow an identical and independent binomial distribution having the parameters 
g{k),pk = cos2^(^. Let rjk = E(k,g{k)) = (^f^) (Pfc)^! " Pfc)^^'^-* = Proh[Xk < 

^ ) • Consider X^^j ~ iid Bernoulli{pk)., where 1 < i < g{k). So, Xj^ = X^^i + . . . + Xj^^gf^j^y 
Thus, var{Xk) = al = var{Xk^i + ...+ Xfc,g(fc)) = 9{k)pk{l - Pk)- 

Since g{k) oo as k ^ oo, the Central Limit Theorem is applicable. Define = ^p^M. 

So r/fc = Proh{Xk < ^^^^) = Prob{Zk < Cfc), where Cfc = Suppose, lim^^oo Cfc = 

C G M. Since convergence in distribution holds in this case, limjfc_»oo % = ^(Oi where <I> is 
the cumulative distribution function of the standard normal variate. Hence, C, = lim^^oo Cfc = 
limfc^oo ^^^^^2_Pfc) 2 ^ limfc^oo ^f?!? = limfe^oo \/ g{k){l - 2pk) as limfc^oo \/pfc(l - Pk) = 

^. Thus, 



c 



limfc^oo 


^gik){- cos 


— limfc_. 


■oo Vgik) sm ^^^^ 


- linifc^ 




— lim^^ 


■<x> 4fc+2 \/9{k). 



So we get, C = when g{k) = t, = —j when g{k) = k"^ and C, = —oo, when g{k) = k^,s > 
2. Hence, lim^^oo % is = <^(0) = 0.5, when g{k) = k, = <I>(-f) > 0.2 when g{k) = k"^ and 
= oo) = when g{k) = k^,s > 2. This gives the proof. ■ 

Theorem 9 Consider Boolean functions on n variables and let N = 2"'. After k iterations, 
the best known classical probabilistic algorithm can distinguish two Boolean functions of weights 
[A^sin^ 2k^i f 1 and [A^cos^ 2fc+i f"! ^^^^ success probability ps = 0.5 when g{k) = k, ps < 0.8 when 
g{k) = k"^ and = 1 when g{k) = k^, s > 2, for large k. 
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k 


The roots of — (top line) and — (bottom line) 


1 


0.250000, 
0.750000, 


2 


0.095492, 0.654508, 
0.345492, 0.904508, 


3 


0.049516, 0.388740, 0.811745, 
0.188255, 0.611260, 0.950484, 


4 


0.030154, 0.250000, 0.586824, 0.883022, 
0.116978, 0.413176, 0.750000, 0.969846, 


5 


0.020254, 0.172570, 0.428843, 0.707708, 0.920627, 
0.079373, 0.292292, 0.571157, 0.827430, 0.979746, 


6 


0.014529, 0.125745, 0.322698, 0.560268, 0.784032, 0.942728, 
0.057272, 0.215968, 0.439732, 0.677302, 0.874255, 0.985471, 


7 


0.010926, 0.095492, 0.250000, 0.447736, 0.654508, 0.834565, 0.956773, 
0.043227, 0.165435, 0.345492, 0.552264, 0.750000, 0.904508, 0.989074, 


8 


0.008513, 0.074891, 0.198683, 0.363169, 0.546134, 0.722869, 0.869504, 0.966236, 
0.033764, 0.130496, 0.277131, 0.453866, 0.636831, 0.801317, 0.925109, 0.991487, 


9 


0.006819, 0.060263, 0.161359, 0.299152, 0.458710, 0.622743, 0.773474, 0.894570, 0.972909, 
0.027091, 0.105430, 0.226526, 0.377257, 0.541290, 0.700848, 0.838641, 0.939737, 0.993181, 


10 


0.005585, 0.049516, 0.133474, 0.250000, 0.388740, 0.537365, 0.682671, 0.811745, 0.913119, 0.977786, 
0.022214, 0.086881, 0.188255, 0.317329, 0.462635, 0.611260, 0.750000, 0.866526, 0.950484, 0.994415, 



Table 1: Roots of = and bk = 0. 

Proof: The proof follows from = 1 - (^^^^(cos^ ^(afc+i) 2(2fc+i) )^^^^~' 

results in Proposition |H1 ■ 
Based on the results of Theorem [7| and Theorem 1^1 it is clear that when the quantum algorithm 
can achieve a success probability almost unity, then the best known classical algorithm can achieve 
a success probability almost 0.5 (from above) after k many steps. The classical algorithm can 
achieve a success probability almost unity only after k'^ many steps for s > 2. Thus the quantum 
algorithm can achieve a quadratic speed up in this case. 

3.1 Distinguishing Boolean functions from two different Sets of Weights 

Let us consider two different sets and Bk, where contains the roots of = and Bk contains 
the roots of hk = 0. From Proposition |31 it is clear that A^ = {sin^ ^^^lyf 1"^ ~ 1, . . . , A;} and 

Bk = {cos^ ^^fiyf 1"^ = 1, . . . , A;}. One can see Tabled for a few examples. 

Given any k, Algorithm |31 can distinguish whether a Boolean function is either from the weights 
yN siv? ^(2T+i)^ f "I from the weights \_N cos^ ^(2™+!)^ f "I ' "^i^h good success probability for m = 
l,...,k. This also helps in solving the question can we distinguish two Boolean functions of 
weights Nwi and Nw2 where wi + W2 7^ 1? As example, one can check the case for /c = 10 in 
Tabled where we may be able to distinguish Boolean functions with weights Nwi and Nw2, for 
wi ~ 0.388740 and 102 ~ 0.317329. However, we leave this for future research as in this paper we 
mainly focus on distinguishing functions having weights Nw and A^(l — w). 

4 Sure Success Weight Decision Algorithm 
4.1 Motivation 

By Theorem 13 we know that the Algorithm is a randomized one and not sure success algorithm. 
In the Grover-like database search algorithms, number of ideas have been proposed to achieve 
sure success. We try to exploit similar strategies here, though we need to make certain subtle 
modifications for this purpose. 
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4.2 Brassard's Sure Success Database Search 

Many algorithms have been proposed for a sure success database search algorithm based on Grover 
search |1 11 llUllTHl 1^ 1131 0]. Meanwhile, we need to reformulate and to generalize the original Grover 
operators, G„, and f/j, for a sure success method as follows. 

Gn = -I\^,){e) = -{I - (1 - e*^)|V^o)(V'o|}. 
Uf = I\solution){4>) = / - (1 - e'-'t')\solution){solution\. 

A sure success database search can find one of solutions exactly by changing two phases, 6 and (j) 
of ■^IV'o>(^) and I\ 

solution) i'P) ■ Can we apply this kind sure success approach in the Grover database 
search into our problem? We try to exploit Brassard method |2] for this. The method is based 
on the following approach. The required minimum number of Grover operation is calculated, 
which assumes to be k. Then, from 1** to {k — 1)*^ operation, the original Grover operation 
is applied. However, for the last A;*^ operation, a slightly different Grover operation is applied 
by controlling two phases, 6 and cf) for and I\soiution){4')- Figures 3(a) and |3(b)1 show 



the Brassard's method in the Hilbert space and Bloch sphere, respectively. Please note that the 
inversion operators, {(^)I\soiuU(m) i'P) in the Hilbert space, correspond to the rotation operators, 

— e**^2 + 2)i?|^^^(— 0) in the Bloch sphere ^2j. As shown in the figures, the last step 
of the Grover operation uses different phases for two inversions in the Hilbert space and for two 
rotations in the Bloch Sphere. 
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4.3 Approach 

At the first sight, it looks that Brassard method can be directly applied for the weight decision 
problem. However, it is a little bit different from the Brassard's case. Brassard method changes 
only the last step because the goal of this approach is to rotate the {k — 1)*^ state to the solution 
state. However, in the weight decision case, we need to satisfy that two initial states for different 
weights should be rotated to the solution and the non-solution state exclusively after operations. 
In other words, if the proposed method rotates the initial state for wN weight case to the solution 
state, the same operation should rotate the initial state for {l — w)N weight case to the non-solution 
state. By this condition, we can make a relation between the initial states for wN and (1 — w)N 
to the solution and non-solution states. Finally, we can decide the weight exactly from measuring 
the final state and evaluating the function with the measured value x. To satisfy this condition, we 
proposed a method which changes the very last two operations. As a result, from 1** to {k — 2)*'* 
operation, we use the {it,tt) phase angles for the and I\soiuUon){4>) operations. However, for 

{k — Vf^ and /c*^ operation, (— ^i, vr) and (—^2, t^) phases should be used, respectively. Algorithm 
m describes the overall idea. 

Algorithm 4 Sure Success Weight Decision Algorithm 

~1. IV'o) = |0)®"|1), i = 0, 

if < Nmin{w, I — w) < Nsin'^^, k is 2, 

otherwise, k satisfies Nsiv? f < Nmin{w, 1 — w) < A^sin^ 2A;+i f • 

2k— 1'^ < P for smaller weight — 2fc+T^' 

^ ~ 2k-l^ ^ P for bigger weight ^ — 2k+l'^ ' 

2. while(i < (/c - 2)) do 
{ 

IV'j+l) = -I\il)o){T^)I\solution){'^)\'4^i) 

i = i + l 

} 

3. \lpk-l) = -I\,Po){-Gl)I\solution){'^)\'4^k-2) 

4. IV'fc) = -I\4,o){-02)I\solution)iTr)\i^k-l) 

5. measure {ipk) in the computational basis, 
let the result be x. 

6-1. if A; is odd and if f{x) = then 
wt{f) = Nmax{w, 1 — w) 
else wt{f) = Nmin{w, 1 — w). 
6-2. if k is even and if f{x) = 1 then 
wt{f) = Ninax{w, 1 — w) 
else wt{f) = Nmin{w, 1 — w). 



4.4 Alternative Final State 

In the Brassard database search algorithm, the final state should be the solution state. However, 
in the weight distinction case, two final states after k operations should be located in the solution 
and non-solution state exclusively. Meanwhile, there are no restrictions of the locations of final 
states except that two final states should be located at solution and non-solution states exclusively. 
Hence the locations of two final states can be alternatively changed with the required number of 
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operations. As a result, we need to think about the final states of the proposed method. The 
previously proposed method, Algorithmic is a special case where the (tt, vr) phase angles are used 
for the /|^(,) (6) and I\ 

solution) i'P) operations. Hence we can infer that the final state for the smaller 
weight case after k operations in the special weight decision algorithm, in Bloch sphere, is 



sm{2k + 


-cos{2k + l)(3 



(1) 



where /3 



2k+l 



vr. Note that all state vectors are represented in the Bloch sphere hereafter. 



Therefore, we can know that the final state for the smaller wight is alternatively changed with the 
required k operations such as 



sin((2/c + l)2^7r) 


-cos((2/c + l)2^7r) 



sin(/c7r) 


— cos(A;7r) 















< 




1 









(2) 



Meanwhile, the (0,0, +1)* and (0,0,-1)* states in Bloch sphere represent solution and non- 
solution states in Hilbert space, respectively. In summary, if the required number of k is odd, 
the final state for the smaller (bigger) weight case should be located in the solution (non-solution) 
state. On the other hand, if the required number of k is even, the final state for the smaller (bigger) 
weight case should be located in the non-solution (solution) state. On the other hand, for the bigger 
weight case, we can easily analyze the alternative final state based on the number of operations 
with the initial angle, vr — /3. This analysis means that we need to find two phase conditions with 
the required number of operations. 



4.5 Modification of Last Two Operations 

Figure |1] explains how we can rotate two initial states to different final states, which should be 
located in the different poles in Bloch sphere as solution, (0, 0, +1)*, and non-solution, (0, 0, —1)*, 
state. Note that Figure |1] shows the case when only two operations are sufficient to decide the exact 
weight. In the figure, the circle and diamond mean two states, which have the smaller and the bigger 
weight, respectively. Our purpose is to find phase conditions, which can rotate two initial states 
to different poles exclusively with the same phase conditions. Hence if the weight is the smaller 
(bigger) one and the method rotates the initial state to the solution state, the method should rotate 
the initial state of the bigger (smaller) weight case to the non-solution state. Note that all figures 
hereafter are viewed from direction in the Bloch sphere for easy understanding. Therefore, 
only X and Z directions and the X — Z plane are shown in all the figures. In the initial step, two 
initial states \ Aq) and |i?o) are (sin/3, 0, — cos /?)* and (sin(7r — /?), 0, — cos(7r— /?))* = (sin/?, 0, cos /?)*, 
respectively. At the first step, two initial states are rotated to \ A\) and \B\) states by using rotation 
R\soiution) (tJ")- Therefore, only the sign of x is changed. In the second step, two last states are rotated 
to \A2) and \B2) states by using rotation R^^^-^^{9i). Meanwhile, \ Ai) state should be rotated to \ A2) 
state, which is the cross point between line Al and line /. The line / is a path, where the point 
l^i) can be rotated by the rotation i?|^jj^(0i). As the same rule, the same rotation operator rotates 
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-2=1^)^ Linel-f 



Figure 4: Last Two Operations 



\Bi) state to I-B2) state, where I-B2) state is the cross point between hne Bl and hne 1 — /. Because 
two last states are rotated with the same phase angle, the value of x is the same, but values of y 
and z have different sign values between them. At the third step, R\soiution){'^) operator is used for 
the last two states. Therefore, the signs of x and y are changed as shown in {A^) and {B^) states. 
At the last step, the two final states are rotated to two different poles exclusively by using the 
same rotation angle, 62. In other words, 1^3) state moves to the —Z pole, non-solution state, and 
|i?3) to the +Z pole, solution state. Finally, if we measure the final state and the measured value, 
then X, is one of solutions (non-solutions) and we can decide that the weight is bigger (smaller) 
one. Meanwhile, the key point of this approach is to find two cross points, IA2) and I-B2) with the 
required number of operations. 



4.6 Correctness 

In the proposed method, we have to change two phases only for the last two operations, not for 
other operations because until {k — 2)*^ operation, there is no cross point such as 1^42) and I-B2) in 
Figure m Therefore, we need to show that until {k — 2)*^ operation, there are no cross point, but 
in the (A; — 1)*'' operation, there are two cross points. For easy explanation, we only consider when 
the required number of operation is odd case and the weight is smaller one. Hence k is odd and 
2k \ < /5 < 2k+i ^' Other cases can be proved with the same approach. 



4.6.1 No Cross Point until {k - 2)*" Operation 



Figure 5(a) shows the last state, IV'fc-2)) when {k—2) operations of -R|^o) {'^)R\soiution) i'^) are applied. 
Note that in the odd k case, \1pk-2) state should be located in the right upper part, i.e., {—X, +Z) 
area. Meanwhile, the line Al, which is the line perpendicular to the axis of \ipi-u,) and meets the 
south pole, is Z = — Xtan/3 — 1. The value of x of the cross point between the line Al and the 
circle is xnne = — sin2/3, and the value of x of \1pk-2), Xk-2, is sin(2(A; — 2) -|- l)/3. Therefore, to 
show that there is no cross point until {k — 2)^^ operation, we need to prove that xune is always 
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larger than Xk-2- 



Theorem 10 (No Cross Point until (fc — 2)*^^ Operation) — sm2/3 > sin{2{k—2) +1) /3 , where 

< /5 < 2^^- 

Proof. From the value of (3, we can get the value of —sin 2/3 as —sm{2^^) < —sin 2/3 < 
-sin(2^). Meanwhile, (fc - 2)7r + < (2A: - 3)/3 < (A; - 2)7r + 2|^. Because A; is odd 

in this case, sin((A; - 2)tt + > sm{2k - 3)/3 > sm{{k - 2)7r + ^^)- Finally, -sin(2^) > 

sin(2A: - 3)/3 > -sm{^^). Therefore, - sin 2/3 > sin(2A; - 3)/3. ■ 

4.6.2 First Cross Point in {k — 1)*^ Operation 

Figure [5(b)| explains why, in the {k — l)*'^ operation, there is the first cross point between the line 
Al and the line /. To prove this, we need to show that the value of x of the cross point, xum, 
between the line Al and the circle is always larger than equal to the value of x of IV'fc-i) state, 
Xk~i- Note that the value of xnm is sin 2/3, and the value of the Xk-i is sin(2(A; — 1) + l)/3. 

Theorem 11 (First Cross Point in {k — 1)*'^ Operation) sin 2/3 > sin(2(fc — 1) + 1)/?, where 

< /5 ^ ra^- 

Proof. From Theorem we can get the value of sin 2/3 as sin(2^^) > sin 2/3 > sm{ 2j^^i )- 
Meanwhile, {k — 1)tt < {2k — l)/3 < (A: — l)7r + -^i+i- Because k is odd in this case, sin(A: — 
1)tt < sm{2k - l)/3 < sin((A; - 1)tt + 2^). Finally, < sin(2A; - 1)^ < sin(2^). Therefore, 
sin2/3 > sin(2A: - l)/3. ■ 

4.7 Phase Conditions 

Figure and [7| show the trace of states for the last two operations when k is even and odd, 
respectively. Note that we only consider the Boolean function with the smaller weight because the 
phase conditions are the same for the bigger weight case. At first, to rotate \1pk-2) state to the first 
cross point |^), 9i should satisfy the following equation 



- sin(2A; - 3)/? 

%>(^i) I 

-cos(2A; - 3)/3 



/ cos(2fc-2)/3-(-l)'= cos^a 
' 2sin/3 

y I • (3) 

y - cos(2A:-2)/3-(-l)'' cos/? 



2COS/3 

As a result, 9i should be chosen as a value to satisfy the following equation 

(-1)'= cos/3 - cos2/3cos(2A; - 2)/3 
^ " sin2/3sin(2/t-2)/3 ' ^ ' 

Meanwhile, the value of y for the state \A) is calculated by y = sin0isin(2A; — 2)/3. With the 
very similar approach we can find a condition for the value of 02, which should satisfy the following 
equation 

/ - cos(2fc-2)/3+(-l)'° cos/3 \ 
2sm/3 

^|t/.>(^2) -y 

\ -cos(2fc-2)/3-(-l)'' cos/3 
\ 2 cos /3 / 




(5) 
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(a) No Cross until {k — 2)*'' Operation 




(b) First Cross in {k - l)"* Operation 



Figure 5: Correctness of Modification of Last Two Operations 



Finally, 62 can be chosen, which satisfies the following equation 

_ (-I)''sin2/?(ysin6l2 - (-l)^'sin/3) 
^ ~ cos/3cos2/3- (-l)*^cos(2fe-2)/? ' 



5 Comparison with Quantum Counting Method 

To argue the efficiency of our algorithm, let us refer to the existing works on quantum counting 
[Uinill]- The existing algorithm exploits the period information of Grover iterations. From this 
period information, one can guess the number of solutions. Hence, as like Shor's factoring algorithm, 
this is the task to find the period of Grover operations using quantum Fourier transform. 



15 




-cos(2/(r-3)^„j^ 

Figure 6: Last Two Operations for Even k 

Because we proposed a way to decide the weight from given special two weights, it is meaningful 
to compare time complexity between our approach and the method based on quantum counting 
with the same promise. 

5.1 Quantum Counting 

Counting problem is to find the number of solutions of a given Boolean function. Meanwhile, 
Grover operation shows some kind period patterns with the iteration numbers. By this analysis, 
we can count the number of solutions with quantum Fourier transform as shown in Algorithm 

Algorithm 5 Quantum Counting [HElll] 



Let 


Cp '■ \m 


) 1^) \m) (g) (Gp)'"!^') 


Let 


Fp : \k) 




1. 


|*o) ^ 


W0 W|0)|0) 


2. 




CfI^o) 


3. 


1-^2) ^ 


^1) after the second register is measured (optional ) 


I 


1^3) ^ 


Fp Olio's) 


5. 


/ ^ measure ^'3) (if f > P/2 then f ^ {P - f)) 


6. 


output: 


N s\T?{f'K/P) (and f if needed) 



5.2 Exploitation and Analysis of Quantum Counting for Weight of Boolean 
Function 

We analyze the quantum counting algorithm for three purposes. First, when a weight w is given, 
how we can check the correctness of the given weight with how much time complexity. Second, 
when two weights are given, how we can decide the real weight with how much time complexity. 
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_ j(e,)«,(;r)|((^.„,._,) 

'cos,S„, +cos(2/t-2),S., ' 
2sin/S„. 
y 

cos^., -cos(2«:-2)yff„^ 
2cos /? 

= ff^(;r)|A) 

-coS;S„_ -cos{2k-2)P,. ' 
2 sin 

~y 

cos/?., -cos(2A'-2)/?,, 
2 cos /?.,^ 



Figure 7: Last Two Operations for Odd k 

Third, when n possible weights are given, how we can find the real weight with how much time 
complexity. 

5.2.1 Check Correctness of a Given Weight w 

If a weight is given as w=s'm'^6, how can we check whether this is correct or incorrect and what 
about time complexity? Let / = P6/tt. If / were an integer, there would be two possibilities: 
either / = (which happens if t = or t = N), in which case l^'s) = |0), or t > 0, in which case 
1^3) = a\f) + b\P — /), where a and b are complex numbers of norm I/V2 HOmi. In other words, 
if we assume the value of P as kir/O, the measured value of / should be k. As a result, we can 
easily check whether the given weight is correct or not by measuring l^'s). If the measured value 
is k, the given weight is correct, otherwise incorrect. In this case, the time complexity is 0{P). 
Meanwhile, because we already know the value of in the initial time, we can find the smallest 
integer value of P as /c^. If is ^, then P is just ka. Therefore, the time complexity of this case 
is 0{ka) when 9 is ^. 

5.2.2 Decide Real Weight w from Two Given Weights wi and W2 

From the previous section, we can know that when wi = sin'^6i{w2 = sin?62) is given, we can find 
the required number P for Wi{w2) and the expected measured value as ki{k2)- However, when we 
want to decide which weight is real one, ki and k2 should be different because they are the clues to 
distinguish. Therefore, we need to find the integer value P, which will be used for both two cases. 
Two values of P are Pi = kiTi/Oi and P2 = /c2vr/02- Because we need to execute the algorithm only 
once, Pi and P2 should be the same. Hence, kiTi/Oi should be equal to A;27r/02 and our job is to find 
suitable ki and /c2- Meanwhile, if 9i = vr/ai and 6*2 = '/r/a2, P should be feiai = ^202- Considering 
the time complexity, we need to find the smallest integer value of P by Pg = LCM{ai,a2), where 
LCM is the least common integer multiplier. Then, ki = Pg/ai and k2 = Ps/a2- Finally, with the 
value of Ps, we evaluate the algorithm, and if the measured value of / is fci, we can say that the 
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weight is wi. If the measured value of / is A;2, we can say that the weight is W2- Time complexity 
of this case is LCM{ai,a2), where 6i is vr/ai and 62 is 7r/a2. 

Now we need to compare our method and the above method based on quantum counting when 
wi + W2 = 1- If wi = sin'^{ 4^+2 ) and 1^2 = sin?{ 4^+2 ) = cos^ ( 4^+2 ) , our method, Algorithm 

k fc+1 k 

and 0] can decide which one is real one in k many steps. On the other hand, the method based 
on quantum counting needs P as 4A; + 2 and if / is k, the real weight is wi and if (/) k + 1, the 
real weight is W2- Prom this analysis, we can know that the method based on quantum counting is 
less efficient than our method because the quantum counting method needs to exploit some period, 
which requires several Grover operations. 



5.2.3 Find Real Weight w from n Possible Weights 

We can extend the previous result to a more general case. When n possible weights are given. Can 
we decide which weight is real one? The most important thing in this problem is to find the smallest 
integer value of Pg = LCM{ai,a2, • • • ,an)- Then, as the same method in the previous section, ki 
should be Pg/di- From this analysis, we can find the real weight from the measured value f as ii f 
is ki then the real weight is Wi. In this case, the time complexity is Ps = LCM{ai,a2, ■ ■ ■ , an)- 



6 Conclusion and Open Problems 

In this work, we have investigated the application of Grover operators to distinguish weight of 
Boolean function when two weights are given. Firstly, when we assume that the weight is the 
number of solutions, we found that Algorithm |31 can find the exact weight with almost certainty. 
Secondly, by exploiting the sure success Grover search method, we found a sure success weight 
decision algorithm. Algorithm |1J with modification of the last two Grover operations with phase 
conditions. Lastly, we have compared the proposed method to the quantum counting algorithm. 
Because the quantum counting algorithm needs period information, which requires more Grover 
operations, our method is more efficient than the quantum counting method. 

On the other hand, until this work, we assume that two weights are given such as wN and 
(1 — w)N. How about other cases such as wi and W2, where wi + W2 ^ 1? Moreover, when three 
or more weights are given, can we find the exact weight with the similar approach? This may be 
attempted using the brief idea presented in Subsection 13.11 
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